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In this paper, we develop a new type of slave particle method which is similar to the slave rotor 
model except that the quantum rotor is substituted by a spin one slave particle. The spin-one slave 
particle itself can be represented in terms of Schwinger bosons/fermions. This approach is more 
conveniently applicable to the strongly correlated Hamiltonians with onsite Hubbard interaction 
and resolves the limitations of using the slave rotor model as well as the Anderson-Zou slave particle 
technique. For instance, the mean-field parameters of the slave spin method do not vanish above 
the Mott transition and this approach is smoothly connected to the non-interacting limit. As an 
example, we study the phase diagram of the Kane-Mele-Hubbard model using our current approach. 
In the absence of the spin-orbit interaction, the Mott transition occurs at Uc — 3ti . Several aspects 
of the slave spin method, its gauge theory and various possible mean-field states associated with 
this approach have been discussed. 
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I. INTRODUCTION 

The theoretical proposal for the topological insulators 
(TI) and their experimental observations have drawn 
a huge interest in the past few years^""'. This new 
state of matter hosts several interesting phenomena. For 
instance, TIs have topologically protected helical edge 
states and exhibit spin quantum Hall effect accordingly. 
Most of theoretical studies in this field have been carried 
out on the weakly interacting systems. However, the ef- 
fects of interaction on the properties of TIs remain largely 
unclear. The intuition is that they are robust against 
weak interactions. However, exotic states of matter or 
spontaneously broken phases may emerge due to strong 
inetraction/corelation. Therefore, it deems necessary to 
study the interacting TI (ITI) in order to have a better 
insight into the physics of TIs. The simplest interac- 
tion that one can consider is the onsite Hubbard interac- 
tion. Recently, several authors have studied the phase di- 
agram of the Kane-Mele-Hubbard (KMH) modeP'® both 
theoretically^"^'* and numerically ^^"^^. 

In this paper, the issue of onsite Hubbard interaction in 
the electronic systems has been addressed. We develop a 
new slave-spin technique that is suitable to study a num- 
ber of famous models in the strongly correlated systems 
e.g. Anderson, Hubbard and the KMH models^^"^^. It 
is worth mentioning that the slave-spin model was first 
introduced by de' Medici et. al. in Ref.^'^ and was stud- 
ied further in Refs.^''''^^. However, our slave-spin model 
is different than the previously known slave-spin model, 
though both approaches share the same physical idea. 
The main difference is that in the previous studies, ev- 
ery electron was attached to its own spin 1/2 slave-spin, 
while in our approach, both spin up and spin down elec- 
trons are attached to the same spin one slave-spin. Ad- 
ditionally, the main point of our present work is to rep- 
resent the slave-spin in terms of Schwinger bosons or 
Schwinger fermions. This is a very crucial feature of our 
model as the spin operators are hard to work with. In 



the previous slave-spin studies, authors did not employ 
this decomposition to deal with the slave spins. 

In the present work, our slave spin approach has been 
applied to the last two models. This approach has two 
important advantages over the two other well-known 
slave particle approaches, the Anderson-Zou^'' and the 
slave rotor models^^'^*. Firstly, it is smoothly connected 
to the non-interacting case and has all the desired prop- 
erties. This was missed in the Anderson-Zou slave par- 
ticle technique at half filling. For example, instead of 
a semi-metal phase, Anderson-Zou slave approach gives 
a superconducting phase for the Hubbard model on the 
honeycomb lattice. Secondly, our model can be easily ap- 
plied above the Mott transition and provides the correct 
behavior in the large U /ti limit. However, the slave rotor 
cannot be easily applied in that limit as the mean-field 
parameters vanishes in the simplest approximation. 

II. METHOD 

Electron carries both spin and electric charge. The 
spin up and spin down electrons carry the same elec- 
tric charge but opposite spins. Similar to the slave rotor 
model, the electron operator in the slave spin model is 
decomposed into a charge creating and a spin creating 
parts. This is a natural way to implement the idea of 
spin-charge separation in the deconfined phases. There- 
fore, we have 

4,. = stfl. (1) 

where Sf creates the charge and //^ the spin of electron. 

The f]^ slave particles that carry the spin of electrons are 
usually referred to as spinous. Since the charge at site i, 
has three possible values 0, — e, — 2e, we assign the charge 
to the Sf — ne{i) — 1 eigenvalues of a spin one object, 
where ne{i) is the total number of electrons at site i. 
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FIG. 1. (Color online) Schematic phase diagram of the slave 
spin approach for the KMH model. t2 is the spin-orbit cou- 
pling constant. There are four phases we can identify in 
the slave spin model: 1- Semi-metal (SM) phase, where both 
charge and spin excitations are gapless. 2- Gapped spin liq- 
uid (SL) phase which according to our theory is going to be 
a chiral spin liquid state. In this phase both charge and spin 
degrees of freedom are gapped. 3- In-plane Neel ordering. In 
this phase we initially obtained a gapless spin liquid where 
charge is gapped but spin remains gapless. The spin liquid 
was coupled to a gapless and compact (7(1) gauge field. After 
including the instanton effect, the spin liquid phase undergoes 
a phase transition into the in-plane (XY) Neel ordering. 4- 
Topological band insulator (TBI). In this phase charge gap 
closes, but spin degree of freedom is gapped in the bulk of the 
system. The spin up and spin down spinous have opposite 
Chern numbers and therefore, it corresponds to an insulator 
with gapless helical edge states. 



Doubly occupied sites correspond to S*^ = -1-1, half filled 
sites carry = and empty site correspond to S*^ = — 1. 
Subsequently, and are the ladder operators for the 
spin one object corresponding to the charge of electrons 
that connects different charge sectors. The physical spin 
of doubly occupied sites as well as empty sites is zero, 
however the spin of 5^ = sector can be either -1-1/2 or 
— 1/2 corresponding to spin up and spin down electrons. 
The following table summarizes the relations 

The electron operator in terms of new slave parti- 
cles needs to be project in order to recover the physical 
Hilbert space. This can be done by implementing the 
following constraint on the physical Hilbert space. 

= nHz)-l = 4A^ + 4/a-l. (2) 

As a sanity check, the above constraints yields four 
sites for the dimension of the local Hilbert space. An- 
other way to confirm this result is to consider the gauge 
freedom in the definition of slave particles. It is clear 
from Eq. [1] that the electron operator is invariant un- 
der the following U(l) gauge transformation 
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TABLE I. The local physical Hilbert space that satisfies the 
Sf = rif (i) — 1 constraint. 



fl^exp{i0,)fl , 5±->exp(T*0,)5±. (3) 

That is the both spinous and carry opposite charges 
under the internal U(l). On the other hand, S'^ {S~) 
carries -1-1 (-1) charge under the spin rotation group along 
z axis {S^ charge), because it changes the eigenvalue of 
operator by -1-1 (-1). Due to the gauge invariance of 
the physical Hilbert states, the total charge of spinous at 
site i has to be equal to the up to a shift which can 
be fixed by acting ^ on one physical state. 

Now, let us use Schwinger slave boson (fermion) rep- 
resentation for the slave spin operators, S^'^. Because 
the Hilbert space of the slave spin has three states, we 
can assign a hardcore boson (fermion) with three flavors 
to each charge state (see Table II). To be more precise, 
6i flavor is assigned to the empty site, 62 flavor to the 
singly occupied and 63 flavor to the doubly occupied 
sites. Therefore, the eigenvalue of the Sf operator is 
the difference between the number of flavor three and 
one Schwinger bosons (fermions). Also, the operator 
acts on the states by changing the flavor of Schwinger 
bosons (fermions). The following equations summarizes 
the scheme 

S! = hlM,, - bUh,.. (4) 

From now on, we only consider Schwinger bosons in 
calculations. The Schwinger boson representation enjoys 
the following internal U(l) gauge degree of freedom (dif- 
ferent from that was discussed earlier). 

^Im exp (m,) 5| ^ , TO = 1,2, 3 (5) 

Therefore they all carry the same charge under the above 
internal U(l) gauge field and is equal to their num- 
ber density. Therefore the total number of Schwinger 
fermions has to be fixed. Another way to obtain this re- 
sult is to consider the dimension of the Hilbert space of 
the charge sector which is three. To achieve that result, 
the number of Schwinger bosons has to be fixed and equal 
to one. Therefore, 

bUK^ + blM + bLh.^^l■ (6) 
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TABLE II. The local physical Hilbert space of the charge 
sector that satisfies the n\ — n\ i +n| j +n\ i — 1 constraint. 



In the path integral formalism, the constraints in equa- 
tions [2] and [6] can be implemented using two La- 
grange multiplier fields and integrating over those auxil- 
iary fields. 

III. ONSITE HUBBARD INTERACTION 

The Hubbard interaction acts locally and adds a cost 
to the creation of doubly occupied sites through Urii^^ni^i 
term where ^ is the number of spin a electron at site i. 
Since rii^^rii^i counts the number of doubly occupied sites, 
i.e. the states with 5f = — 1 = 1, it is equal to the 
number of 63, i Schwinger bosons. Therefore n^.^rii^^ — 

&3 j^s^i- So we have 

n^,tn^,i = 4 ,,63,i (7) 
i i 

On the other hand, due to the Sf = n{ — 1 constraint, 
at half filling we have 

E = E (^3/3,. - &U&m) = E - 1) = o-(8) 

i i i 

So we can rewrite the Hubbard interaction as follows 

^u = ^E(^L^3,, + 6lA»). (9) 



IV. GAUGE THEORY OF THE SLAVE SPIN 
MODEL 

As we discussed before, there two U(l) internal gauge 
degrees of freedom. U{l)s represents the gauge the gauge 
freedom in the definition of electron operators and U{l)c 
that represents the gauge freedom in the Schwinger bo- 
son representation of spin slaves. Table III summarizes 
the transformation properties and the charge of different 
slave particles under U{l)c and U{l)s- 

Due to the compact nature of U{l)s and U{l)c gauge 
fields, both fields have to be gapped to obtain a stable and 
physical mean-field results. Otherwise, instanton prolif- 
eration destabilize the mean-field state. Gauge fields can 
be gapped out in two ways: 1. through Anderson-Higgs 
mechanism, 2. by the presence of Chern-Simons action 
in the low energy physics. In the following, we comment 
on the possible outcomes of the slave spin model. 
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TABLE III. The charge of the different slave operators un- 
der electromagnetic U{1)em, U{1)s and (7(1) c gauge fields. 
The transformation of charge q operators after 6i gauge trans- 
formation is Og{i) — >■ exp {i9i) Oq{i). The transformation 
of bl ^ quasiparticles under ?7(1)3 have been chosen so that 
St — exp (=piSi) and Si — > S'f. The electromagnetic 
charge of spinous —Qs can be assumed to be any number. 
This ambiguity will be resolved if we appropriately take the 
U{l)s gauge field into consideration which leads to the gener- 
alized loffe-Larkin formula. For more details see Refs.2^ and^° 



A. U{l)s gapped, Uc gapped 

In this case, gauge field is gapped and instanton oper- 
ator is irrelevant. Therefore, the mean- field calculations 
can be trusted. Four possibilities may happen: 

1- Both spin and charge degrees of freedom are 
gapless. — The spin gap can be attributed to the gap in 
the spectrum of spinous {fi^a particles). The charge gap 
on the other hand is physically related to the doublon 
(doubly occupied site) or holon (empty site) creation. 
Hence, the charge gap can be attributed to the gap in 
the spectrum of 63.^ and bi^i Schwinger bosons (note that 
62,1 represent the singly occupied sites). When instan- 
tons are absent and both spin and charge degrees of free- 
dom are gapless, i.e. there is no gap in the spectrum 
of spinons and 63^^ or bi^i slave particles, the mean-field 
state is nothing but a conducting phase and is a Fermi 
liquid (FL) system. 

2- Spin is gapped while charge is gapless. — This 
phase corresponds to the band insulator, where band- 
structure calculations for electrons predict a filled band 
that is separated from the conduction band by a nonzero 
gap. Although the above statement is valid for the bulk 
of the system, there may be gapless edge states present. 
Two famous examples are the non-interacting Kane-Mele 
and Haldane models, where they host helical and chiral 
edge states respectively, though their bulk are insulating. 

3- Spin remains gapless while charge is 
gapped. — This phase corresponds to the Mott insulat- 
ing phase, where band-structure calculations for electrons 
predict a metallic phase, while taking correlations into ac- 
count changes the behavior of the system. In Mott insu- 
lator, charge cannot be excited due to the nonzero charge 
gap while spin degree of freedom may or may not be gap- 
less. If there is no magnetic ordering in the ground-state 
of a Mott insulator, the system is called a spin liquid. In 
this case, we obtain a gapless spin liquid. 
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4- Both spin and charge are gapped. — This phase 
simply corresponds to the gapped spin hquid provided 
that there is no magnetic ordering in the ground-state. 

B. U{l)s gapped, U{l)c gapless 

In this case, U{l)c gauge field is gapless and instanton 
operator that adds a quantum of flux of the correspond- 
ing field through quantum tunneling is most likely a rel- 
evant perturbation. Under instanton proliferation, the 
mean-field state becomes unstable and we cannot trust 
our results. Instanton derives the system into the confine- 
ment phase where only U{l)c gauge invariant operators 
can be measured and any other operator vanishes. 

C. U{l)s gapless, U{l)c gapped 

In this case, U{l)s gauge field is gapless and again in- 
stantons will proliferate and the mean-field result cannot 
be trusted. 

D. U{l)s gapped, U{l)c gapless 

In this case, both internal gauge fields have gapless 
modes that results in instanton proliferation. The mean- 
field state is again unstable and slave particle is not a 
good approach in this case. 

V. PHASE DIAGRAM OF THE 
KANE-MELE-HUBBARD MODEL 

The Kane-Mele-Hubbard model on the honeycomb lat- 
tice is described as 

(10) 

where ii, U, and t2 are the nearest neighbor hop- 
ping energy, the strength of the on-site repulsion, and 
the second-neighbor spin-orbit coupling strength, respec- 
tively. Here Cia (cl^) annihilates (creates) an electron 
with spin a on site i. is introduced so as to obtain a 
nonzero fiux turning around any triangular path and is 

defined as Vi = .'b^'^J . .z where d, and are two short- 

est vectors that connect sites i and j i.e. di + dj = Ri — Rj 
(see Fig. [2]). 

Now we can employ the slave spin model to rewrite the 
KMH model in terms of spinous and Schwinger fermions. 
To do so, note that 

cI,(TCj,(T = S^Sj fl^crfj.iy- (11) 

where Sf {S~ = S^^ ^) is defined in terms of Schwinger 
bosons in Eq. [4]. Now we can decouple spinous from 




FIG. 2. (Color online) Kane-Mele model for the topological 
insulators. Arrows denote the phase of the hopping to the 
NNN for the spin up electrons. The phase for the spin downs 
is opposite to that of spin up electrons. Assigning a phase ^o- 
to the hopping of spin a electrons induces a 3(j>cr magnetic fiux 
for spin a electrons moving around the depicted triangles. It 
is this nonzero spin dependent fiux that causes spin quantum 
Hall effect for non-interacting electrons. 



slave spins using the following Hartree-Fock approxima- 
tion that can be justified through Hubbard-Stratonovic 
transformation 

4,o-Cj,cr —{S^Sj ) fj^crfj.a- + ^t^j (^fi,afi-o^ 

-{StS-){flj,,.). (12) 

Now let us define the Xa [hj) mean-field parameter in 
the following way 

xAh3) = {flj,.a). (13) 

For the nearest neighbor, we assume a uniform and spin 
independent xi (i, j) as follows 

xi(«,j)-Xi/2. (14) 

For the second nearest neighbor, except in the chiral spin 
liquid phase that will be discussed later, we assume that 
it takes the following value 

xl{i,i) =iVi,j(JX2l'2- (15) 
Using the definition of the slave spin operators we have 

sts- = (blM, + hhK) (4,,&3,, + fol.,62,,) .(16) 

In the Schwinger boson approach, it is more convenient 
to assume ih^ra,i^n,j^ = where 11,171 — 1,2,3. Conse- 
quently, 

StSj ^ A22 blfil ^ + Ai3 (*, jT blfil, 
+i?.c.-2i?e(A22(«,j)A*3(*,j)), (17) 

where 

Ai3(«,j) = (fel,*^3,,), (18) 
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and 



A22 ^ (&2.»&2j) 



(19) 



As a result, 

(5+57) = 2i?e(A22(*,j)AJ3(*,j)). (20) 

We assume A22 («, j) and A13 are real and indepen- 
dent of direction and only depend on the distance of two 
neighbors. The Hubbard term can also be written as in 
Eq. [9] . The microscopic constraints on the Hilbert space 
of slave particles in Eqs. [2] and [6] can also be imple- 
mented at the mean-field level using Lagrange multipli- 
ers Ai and A2. It is straightforward that at half filling, 
Ai = naturally satisfies the constraint in average. Us- 
ing all these procedures, we finally obtain the following 
form for the effective mean-field Hamiltonian 



At half filling, all the negative energy levels are occu- 
pied and the gap in the spectrum is proportional to 1^21- 
Although the bulk is insulating, the edge of the system 
can be shown through various techniques e.g. Jackiw- 
Rebbi solitons''^ host gapless spinon excitations, even in 
the disordered system. The number of gapless modes for 
each spin equals to the |Cct|, where is the first Chern 
number of the band structure. Using the continuum limit 
of the above Hamiltonian that yields two gapped Dirac 
cones can be used to compute the Chern index which is 
going to be = — = C and |C| = 1. 

Now let us comment on the effective Hamiltonian for 
Schwinger bosons which is 



U 
2" 



^eff ~ + Hb + He 



(21) 



A 



where 
Hf 

(22) 
where 



<i,j> 

-2hxi V A 



22 \^''3,i''l,: 



<i,j> 



( fk,A.a fk.B,a ) ( 1 



t*2<^Ck tlrik 
V*k ~i*2(^(.k 




I3 (4/L +^2j^2.) 
2t2X2 ^l2(biA.j+b^,,K 



<<i,j» 
<<i,j» 



A^3(4/2,,+^2,,^2,« 



(28) 



t\ — 2tiA22Ai[3 
t*2 — 2t2^^2^1Z- 



(23) 
(24) 



Since 63, and 61 Schwinger bosons are paired and are 
decoupled from 62 Schwinger bosons at the mean-field 
level, we have 



Hh — Hi. 



H, 



b J 



and the following structure factors have been used 
rjk ~ cxp {—iky) + 2 cos |^\/3A;2;/2^ exp (iky/2) (25) 
Ck = 2sin [^/ik^/2}j [cos {iky/2) ~ cos [V^k,^ /2y^{2Q) 

Thus, the energy spectrum of the spinous are 



in which 



We have used the following definitions 
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^ b 



'l,A.k 



A. 



,k^"3,A 



'3,B-k 



(29) 



(30) 



(31) 



Ef, = ±^/it*aa)^ + \qrjk\' 



(27) 



and 



Ml' 



( U/2-\2 -t2X2^l2ik -hxiAl2Vk\ 

U/2-X2 -hxi^l^Vl ~t2X2^Uk 

-t2X2^l2^k -hxi^hik C//2-A2 

X-hxi^hvl -t2X2^l2S,k C//2-A2 ) 



(32) 



where 



= 2 cos {Vik^/2^ cos {iky/2) + cos [Vik^^ . (33) 
The energy excitation of 63 — bi branch can be easily computed using the Bogoliubov transformation and is given by 



E. 



13, 



{U/2 - \2f ± 2tit2XlX2^h2^l2 \rik\ a - \hXl^\2Vk\ - {t2X2^l2is 



(34) 
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For the 62 branch we have the following pairing Hamiltonian 



22, fe 



Mr = - 



13, fc 

= (4,A,fc' 4,S,fc' b2,A,^k, b2,B-k 

/ A2 2i2X2A?3,fcCfe 2tixi^\^^k\ 

A2 2tiXiA}37;*. 2i2X2A23a 

2t2X2A?3efc 2tiXiAi3r/fc A2 

V2tiXiAi3ry* 2t2X2A23efc X2 J 



(35) 



The energy spectrum of the above Hamiltonian can also be computed through the Bogoliubov transformation and is 
given by following expression 



E22.k = V^l ± 8^1^2X1X2 AI3A23 |77fc|a - \2hXl/^\M - {2t2X2^l^^s,k) 



(36) 



There is also a classical contribution to the energy 
which can be calculated by adding c-numbers in the pro- 
cedure 



12iixiAi2Ai3 + 24t2X2A22A23 + -A2 



U 



(37) 



where Ns is the number of sites. To obtain the varia- 
tional mean-field parameters, we need to compute the 
total ground-state energy and minimize it with respect 
to its parameters. It is given by the following relation 



E- 



E. 



22, fe 



E. 



22, fe 
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VI. POSSIBLE PHASES IN THE 
KANE-MELE-HUBBARD MODEL 

To interpret the results of the slave spin approach to 
the KMH model, we would like to make several comments 

I-Single condensation versus pair condensation 
of bosons. Boson can condense both as local objects 
(singly) or form pairs and condense as extended objects. 
Usually start from pair condensation scenario and solve 
the Hamiltonian. When the energy excitation gap closes, 
it is easy to argue that single condensation scenario is 
energetically more favorable* . In the slave spin approach 
to the Kane-Mele-Hubbard, our numerical study shows 
that the energy spectrum for the 62 Schwinger bosons 
is gapless in the whole phase diagram and therefore we 
have 



2,j 



^0 



(39) 



II- When do internal gauge fields are gapped? 

To gap out a gauge field, we have two options. 1- con- 



dense a charged operator. Any operator Oq{x) that car- 
ries charge q under that gauge field, transforms after f{x) 
gauge transformation in the following way 



dq{x) exp {if{x)) dq{x). 



(40) 



If the system is gauge invariant, the expectation value 
of all physical quantities has to be gauge independent as 
well, i.e. 



6,(a;)) = (exp (j/(x))6,(x) 



(41) 



The above equation has two solutions: q = or 
Oq{x)^ = 0. Now what if we find a charged opera- 
tor whose expectation value is nonzero? It simply means 
that the gauge symmetry is broken. This is the reason 
why Cooper pair formation leads to the Meissner effect 
and as a result to the superconductivity. When gauge 
symmetry breaks, gauge particles acquire mass through 
rin?a^a^ term in their action. It is easy to show (e.g. by 
using the Landau-Ginzberg theory of phase transition) 



that 



2- Through Chcrn-Simons term. When the ground- 
state breaks the time reversal and exhibits quantum Hall 
effect in response to the gauge field that we are arguing 
about, the low energy action for that gauge field is given 
by the following Chern-Simons action 



G.F. 



An ^ 2e2 



{e^^Pd^axf (42) 



It is straightforward to see (e.g. by solving the classical 
equations of motion for the gauge field) that the gauge 
particles acquire a nonzero mass proportional to |C|. 

As we mentioned above, our numerical study for the 
KMH shows that always < >^ 0. Since b\ i carries 
nonzero charge under the internal [/(l)c, its condensation 
leads to the gauge symmetry breaking of the C/(l)c gauge 
field. We only need to be concerned about the U{\)s 
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gauge field. Our numerical and theoretical observation 
at half filling, and when i2 = is as follows: 

U < Uc.i — 3ti . — The energy excitation for the 
63 and 61 vanishes below ?7c,i and they undergo Bose- 
Einstein condensation. Therefore we have 



t2 - 0, U < Url 



^0. (43) 



The above equation simply means that the U{l)s gauge 
symmetry breaks as well below Uc,i- At t2 = 0, spinous 
are gapless as well, and altogether we conclude the system 
is in the semi-metal phase. 

Uc,i < U < Uc.2 ■ &3.1 Schwinger bosons are 
gapped, however spinons form chiral state with 
nonzero total Chern number. — This result needs 
more explanation and is beyond a simple meanfield re- 
sult. When the charge gap is nonzero (63 and bi 
Schwinger bosons form bound-state), the charge degree 
of freedom freezes and can be safely integrated out to find 
effective action for the spinon degree of freedom. In our 
case, it corresponds to the loop corrections to the spinon- 
spinon interaction. Such approach has been studied by 
one of us in Ref.® in detail. When t2 = 0, the effective 
Hamiltonian for the spinons is given by J1—J2 Heisenberg 
model. In Ref.^ authors have considered the more general 
case where spin-orbit interaction is taken into consider- 
ation as well and we obtained the so called Kane-Mele- 
Heisenberg model. They have obtained a region above 
Mott transition where hosts spinons with gapped chiral 



ground-sate and nonzero Hall response (nonzero Chern 
number). The nonzero Chern number of spinons in re- 
sponse to the U{l)s gauge field opens up a nonzero gap 
in the spectrum of gauge particles and breaks the gauge 
symmetry. Therefore in this phase, both internal gauge 
fields are broken and we do not need to worry about in- 
stanton effect accordingly. Since, both charge and spin 
gaps are nonzero and spinons form a chiral state, we 
call this state gapped chiral spin liquid which extends to 
nonzero t2 as well. 

U > Uc,2 ■ — In this case, charge gap is again nonzero, 
while J2 is very small and is not enough to form chiral 
state for spinons. Therefore, spinons are gapless in this 
case. Since non of the Anderson-Higgs mechanism, nor 
the Chern-Simons action break the U{l)s gauge field, in- 
stantons will proliferate. Instanton effect results in spon- 
taneously broken phases such as Neel order or valence 
solid bond (VBS) orders®'^^. 



VII. CONCLUSION 

In this paper, we introduced slave-spin model which is 
similar to the slave rotor model. It was discussed that 
using this technique is more conveniently applicable as it 
is applicable in the whole phase diagram and is smoothly 
connected to the non-interacting electron systems. The 
gauge theory of this model has been presented. It was ar- 
gued that this approach explains the existence of gapped 
spin liquid in the phase diagram in a narrow region for 
moderate values of Uc/ti and for small values of ^2- 
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